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Abstract 

This paper studies the second-order traffic flow models from t he aspect of data enabled 
computations. We consider the phase transition mode l (PTM) ( Colombo! . 120021) as well 
as its extension with a Siebel-Mauser type source term ( Siebel and Mauserl l2006ah . under 
various assumptions of traffic conditions. We derive the estimation of first-order and 
second-order traffic quantities based on the vehicle t rajectory data. The dataset we utilize 
is the Next Generation SIMulation (jNGSIMI . I2006T ) which contains high quality vehicle 
trajectories on a segment of 1-80, California. The NGSIM data records the location of 
each passing vehicle at a high precision of every 0.1 second. Our goal is to perform 
traffic estimation using different second-order fluid models, and analyze how the quality 
of estimation deteriorates with less frequent sampling of vehicle locations. The findings 
of this paper will help evaluate the effectiveness of current mobile sensors such as GPS 
or smart phone - usually with a sampling period of 3 seconds - in reconstructing traffic 
states. 



1 Introduction 



1.1 Aw-Rascle-Zhang model and the Siebel-Mauser type source term 

The Aw-Ras cle-Zhang model was independently developed by Aw and Rascle ( 2000l ) and 



Zhand (|2002h . The model is established as a two-by-two system of conservation laws. The 
first equation is based on the conservation of mass: 



pt + (vp) x 







(1.1) 



for now and sequel, subscripts t and x denotes partial derivative with respect to time and 
locatio n, respectively. Equation (11-11) is similar to t he classical Lighthill-Whitham-Richards 
model (jLighthill and Whithaml . Il955l : lRichardd . fl95i b but without assuming that the average 
velocity v depends only on the average density p. The second equation of the model can be 
derived from a car following model. Taking the reference frame of a driver, and denoting x{t) 
its position and v(t) its velocity, we write the equation: 

d r , u(p(t,x(t)) - v(t) 



dt 



[v{t)-u(p(t,x(t))} 



1.2 
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where the equilibrium velocity u depends only on the average density p. T G M+ is a relaxation 
time. In Eulerian coordinates we write: 



(v -u) t + v(v- u) 3 



u 



(1.3) 



and using u t + vu a 



-pu'v x we derive: 



v t + (v + pu'(p)) 3 



u(p) 



;i-4) 



which is the second equation. 

The second equation is then further mani pulated algebrai c ally t o reach the form of an 
homogeneous systems o f conservation laws, seelAw and Rascle (120001 ) for more discussion. 

Siebel and Mauser (ISiebel and Mauserl . l2006al ) proposed a variation of the Aw-Rascle- 
Zhang model by taking into account the reaction time of drivers. The construction is the 
same as above but gives rise to a second equation as: 



v t + (v + pu'(p)) x 



u(p) 



;i.5) 



where r is a reaction time. Typically r (in second) varies within [ 0.5, 1], One should c omput e 
T as the time at which the velocity is essentially at equilibrium. In lSiebel and Mauserl (|2006al ). 
T was chosen to be 2/3 second. 

Finally T — r is modeled as a factor which can be positive (for very small or very high 
densities) or negative (for intermediate densities). One has to notice that a negative factor 
gives rise to stable traffic, while a positive ones produces instabilities. 



1.2 Phase transition model 

The h yperbolic phase transition model (PTM) for traffic flow was introduced by Colombol 
(|2002h . ~.t was motivated by the empirical observation that when the vehicle density exceeds 
certain value, the values of the density-flow pair are scattered in a two-dimensional region, 
instead of forming certain one-to-one relation as is assumed by the classical LWR model and 
the fundamental diagram. 

The PTM consists of two phases: the free flow phase and the congested phase. In the free 
flow phase the dynamic is governed by the LWR model 



Pt + [p ■ v] 3 
v = v(p) 







;i.6) 



The congested phase is governed by the following system. 

' p t + [p-v] x =0 
Qt + [(q ~ q*) ■ v] x = ( 
v = v(p, q) 



1.7) 



where q is identified as the momentum, q* is a given parameter, v is assumed to be a known 
function of both the density p and the momentum q. One choice for such a function is 



v(p, q) 



1 



P 



Pjo 



(1. 
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where Pj a m is the jam density. Another choice is 

v(p, q) = A(p jam - p) + B{q- q*)(p jam - p) 



;i.9) 



for some given parameters A and B. 

One has a lot of freed om in choosing the second order model in the congested phase. 
For instance, Goatin (|2006l ) proposed a phase-transition model which employs the Aw-Rascle- 
Zhang equations in the congested phase. In particular, we can consider the PTM with modified 
congested velocity by adding a relaxation term as in the Siebel-Mauser model. 



1.3 The Next Generation SIMulation (NGSIM) 

Initiated by the United States Department of Transportation (US DOT) F ederal Highway 
Administration (FHWA) in the early 2000's, the Next Generation SIMulation (jNGSIMl . 120061 ) 
program collected high-quality primary traffic and trajectory data intended to support the 
research and testing of the traffic models and open behavioral algorithms. 

We take the 1-80 data set collected and processed on a segment of Interstate 80 located in 
Emeryville, California on Apr. 13, 2005. A total of 45 minutes of data are available, segmented 
into three 15 minute periods. The data set contains vehicle trajectory, instantaneous velocity 
and acceleration recorded at a high precision of every 0.1 second. By neglecting sensing and 
processing error, we treat the raw data set as the ground-truth traffic state. 



2 Traffic Data Measurements 

Using onboard devices such as GPS or smart phone, it is possible to measure the position 
of a car every St seconds. Here 5t is related to the device's characteristics such as desired 
precision and transmission capacity. Let us denote x(t) the position and v(t) the velocity 
at time t. Assume we measure the position at three consecutive shots t\,t2 and t 3 with 
t2 — t\ = t 3 — t2 = 5t (say around 3 seconds for GPS). From these measurements we can 
deduce the approximate velocities in the intervals [£1,^2] and [£2^3] as: 

x(t 2 ) -x(h) x(t 3 ) - x(t 2 ) , . 

«i,a = If , ^2,3 = ^ • (2.10) 

The velocity at time £2 is approximated as 

V(h) = — (2.11) 

One also gets 

D r \ ^23 — ^12 x(tr,) — 2x(to) + x(ti ) 

where D/Dt = d/dt + v ■ d/dx denotes the material derivative in eulerian coordinates corre- 
sponding to the acceleration of the car in Lagrangian ones. 

Another important quantity to estimate is the spatial variation of velocity in Eulerian 
coordinates. For shorter notation, we denote Xi = x(ti). Assuming a mild variation in time of 
the eulerian velocity v(t,x), we write: 

( x 2 + xA ( x 3 + x 2 \ 

V t 2 , ~ Ul 2 , V t 2 , ~ t>2,3, 
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from which we get, setting Sx = X3 + X2 — x 2+ x i ■ 
in other words: 

d~x V [ h > 4 J ~ ~^T- = J t X3 - Xl (2 ' 13) 

Clearly such approximation is acceptable as long as the variation between ux,2 an d i>2,3 is not 
too large. 



3 Traffic Quantities and Model Fitting from Mobile Data 

Depending on the choice of the model and the level of approximation we can use the obtained 
traffic data in many ways to fit an evolutive model (with Kalman filter or other methods). 
Let us describe in this section some choices which can be used depending on traffic data type, 
level of accuracy and computational constraints. 



3.1 Phase transition model 

We depart from the basic phase transition model introduced in Section 11.21 First note that 
we should determine two critical velocities Vf > V c and use either LWR model (i.e. the model 
corresponding to the free flow phase) if v > Vf, or the second order model if v < V c . In the 
first case the velocity is sufficient for data fitting. In this paper, we will focus on the second 
case. 

Remark 3.1. If one wants to use a Newell-Daganzo type fundamental diagram (with constant 
velocity in the free flow phase), then the equation for v is trivial, but the density can not be 
reconstructed. 

3.1.1 Phase transition model with source and strongly stable traffic 

We fix the phase transition model (|Colombd . l2002T ) with velocity in the congested phase given 



by 

V(p, q) = MPjam ~ P) + B(q - q*)(pjam - p) (3.14) 

In addition, we add a Siebel-Mauser type source term. Thus the second equation is given by: 

qt + (v(q-q*)) x = i 4 F Z ^- (3-15) 

1 — T 

If the traffic is stronly stable, dur ing the data measurement we may assume pt, Px,1x ~ with 
only qt non- vanishing. Following ISiebel and Mauser (|2006ah . we assume T ~ 2/3 second and 



r ~ 1 second. Then T — r ~ — l/3sec. Notice that by assumption, f|3. 15j) reduces to 

Qt = 7^ (3.16) 
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Then we can derive 

— v = d t v(p(t, x), q(t, x))+v- d x v(p(t, x), q(t, x)) (3.17) 

= v p ■ p t + v q ■ q t + v ■ (vp ■ p x + v q ■ q x ) (3.18) 
= v q -q t (3.19) 

= B(p jam - p) ■ £^-£ (3.20) 
= - 3B(p jam - p)(q - q*) = - 3 (v - A(p jam - p)) (3.21) 

From the strong stability we get: 

{q-q*)t = f^ = -3(<Z-g*) (3.22) 
Taking into account only the measurements of v and Dv/Dt, we deduce from (|3.2ip that 

1 ( lDv\ 

(P jam -P) = 2 [v+--) (3.23) 
1 A — 

('-"')= -sbU^E (3 - 24) 

Using the same notation as in Section [21 the quantity Dv / Dt is computed approximately at 
time t2', while to get a value at t-i of v with good approximation one should set vfe) ~ v ^^+ v ^,3 _ 

3.1.2 Phase transition model with source and less strongly stable traffic 

In this case we are relying on less strong stability of traffic conditions. In other words, we 
no longer assume that p x vanishes, while still neglecting pt and q x . Then the equation (I3.15|) 
becomes 



qt + v x (q-q*) = ^—^ (3.25) 

1 — T 



We can now write 



^ = v t + vv x ~ v q q t + vv x = v q f j; - (q- q*)v x ^j + vv x (3.26) 



By introducing the variables p = pj am — P, Q = Q — Q* , we obtain 

v = A{p jam - p) + B(q-q*)(p jam - p) = p(A + Bq) (3.27) 
and deduce from (|3.27p that 

Dv 



D/ VV X = V q ~ V Vx ) = B P ~ ^ ' ^' 28 ^ 

From (I3.27p . (13.280 we immediately get the expressions for p and q in terms of v, v x and j^. 



A {T-r){%-vv x ) 



(3-30) 



Dt 

v x are given by (12,121) and (|2.13p respectively. To approximate v, notice that if one 



chooses (|2.10p . the quantity y£ — vv x will always vanish, making q zero. Thus to avoid such 
trivial solution, one should approximate v by either 

x{t 2 )-x{ti) x(t 3 )-x(t 2 ) 
v{t 2 ) ~ v h2 = ft , or v(t 2 ) ~ v 2}3 - 



St K " ^° St 

3.1.3 Phase transition model without source term 

In this case, we consider the PTM without the Siebel-Mauser type source term: 

Pt + (p- v) x = 

qt + {{q-q*)-v) x = (3.31) 
v{p, q) = A(p jam - p) + B(q - q*)(p jam - p) 



We start with the identity 

Dv 
~~D~t " 

Using (|3.3ip . we deduce 



vt + v ■v x = v p pt + v q q t + v{v p p x + v q q x ). 



Dv 

— = v p (-p x v - pv x ) + v q (-q x v - (q - q*)v x ) +v(v p p x + v q q x ) 
= -v x (v p p + v q (q-q*)). 

If we consider the measurements of v, Dv/Dt and also v x as indicated in Section [2] then we 
can consider the two relations: 

v = A(p jam - p) + B(q - q*){pjam ~ p) (3.32) 

^ = v x (Ap + B(q - q*)(2p - p jam )). (3.33) 

Recall the variables p = Pj am — P and q = q — q*. Combining (|3.32|) and (|3.33p and solving 
for p we get: 

P =B-q{ V + - x m- Ap H ^ 
So q and v + ~ Apj am always have the same sign. Substituting (|3,34p into (|3.32[) . we get 

B 2 p jam f + B^2Ap jam -^^-2v^q-A^v + ^^-Ap jam SS j =0 (3.35) 

The discriminant of the above quadratic equation in the variable q is 

.„ 2 ( 1 I Dv . \ 2 , . „ 2 / 1 Dv . 
A = AB Z I v + - — — - Ap jam \ +4AB z p jam [v + —— - A Pjam 

2 
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In order for any meaningful real root of (|3.35p to exist, a necessary condition is that 

1 1 Dv 



v + 



2v x Dt 



> \[Ao~ (3.37) 



In the case where the strict inequality of (|3.37p holds, two distinct roots q\ and tj2 exist. We 
distinguish between two cases: 

• If v + — Apj am > 0, then (|3.35|) has one positive root and one negative root. By 
(|3.34|) . one should choose the positive root since p must be non-negative. 

• If v + 7^757 — Apj am < 0, then (|3.35[) has too roots with the same sign. According to 
(|3.34p . both roots should be negative. This requires that 



1 Dv 

2Ap jam — — — 2v > 

v„ Dt 



A sufficient condition for the above to hold is that -7- -7^7 > 0. We will discuss this 
condition more in the numerical study at Section SJ 

3.2 Estimation of mission and fuel consumption 



Following iBarth et al. (|2007h . the emission and fuel consumption can be modeled via the 



power 

demand function. By knowing the vehicle's attributes such as mass, cross-sectional area and 
rolling resistance coefficient, and the prescribed acceleration and velocity on a particular grade, 
the total tractive power requirements (in kilowatts) placed on a vehicle is given in simplified 
form: 

ptractive = J£_ . . / + . ^ m + ? M . . q + P . 2 . A . q \ . ^ ^ 

1000 y 1 v 2 ; 1000 v ' 

M = vehicle mass(1200 kg) 

v = vehicle velocity (m/s) 

a = vehicle acceleration (m/s 2 ) 

g = gravitational constant(9. 81 m/s 2 ) 

9 = road grade angle 

C r = rolling resistance coefficient(0.005) 

p = mass density of air(l. 225 kg/m 3 , depending on temperature and altitude) 

A c = cross-sectional area(2.6m 2 ) 

C a = aerodynamic drag coefficient (0.3) 



The values for certain parameters above are taken from Chan ( 19971 ). For simplicity, we assume 



all vehicles under study have the same attributes, e.g. mass, cross-sectional area, etc. The 
vehicle emission or fuel consumption can be modeled simply as an affine function of p tractlve ; 

Emission = (3 + aP tractive (3.39) 
Fuel Consumption = b + aP tractive (3.40) 

for some constants a, a > 0, f3, b. 
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4 Numerical Study 



The goal of this section is to numerically verify whether data collected via mobile devices such 
as GPS or smart phone are sufficient to reconstruct first- and second-order traffic quantities. 
To this end, we will utilize the GNSIM vehicle trajectory data in the following way. 

First, we reconstruct the GPS data from the raw data, which records vehicle location every 
0.1 second, by taking samples at every N data points. Here N is a positive integer. Since 
the sampling period of mobile devices is usually 3 seconds, a typical value of N is 30. For 
comparison purpose, we consider other values of N as well. 

Second, we use the reconstructed GPS data from previous step to approximate basic 
traffic quantities discussed in Section [2l More specifically, let ij, i = 1, . . . , m be the time 
instances at which the location is reported by the GPS, where ti+i — t{ = 5t ~ 3 seconds. 
Let x(ti), i = l,...,m be the corresponding location of the vehicle. The estimates of the 

Dv 

Dt ■ 



quantities of interest v, jS, v x read 



V( tl ) ~ Sfc+l)-^-!), or v{U) „ *(*<+!)- (441) 

26t ot 

m {u) && (442) 

d_ f x(t i+1 ) + 2x(tj) + x(t i - 1 ) \ 2_ x{t i+1 ) - 2x(tj) + x{tj-i) 

dx v y u 4 ) ~ st' X {t i+l ) - xit^) [AA6) 

Third, using the estimated values of v, j£ and v x in the previous step, we perform the 
estimation of p = pj am — p (first-order) and q = q — q* (second-order) as in Section EH for 
the following cases. For shorter notations, we denote by v, j£ and v x the right hand sides of 
(|QT|) . (g32D and ([Qg]) . respectively. 

• Phase transition model with source and strongly stable traffic: 

p{ti,x(U)) = LL+l^j (4.44) 

3 Dt 

Here in (|4.45p . the velocity v is estimated using either expression in (|4.4ip . in order to 
avoid vanishing denominator v + | . 

• Phase transition model with source and less strongly stable traffic: 

v + lDv 

xM > ~ v m^7, <446) 

*(*<)) = - UffiSt iiA7 i 

u 3 Dt 

Notice that if we approximate v with the first expression of (|4.4ip and approximate 

Dv 
Dt ■ 



§|, v x with (14421) . (14431) respectively, then 



Dv x(t i+ i) - 2x(ti) +x(U-i) 

Wt ~ = 

Dt x St 2 

_ x(t i+ i) - z(ti-i) _ 2_ _ x(t i+1 ) - 2x(tj) + x(tj-i) = 

26t ' 5t' x(ti+i) - x(ti-i) - ■ y- ) 
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Thus, to avoid such trivial solution for q, we will use the second expression of (|4.41|) for 
estimating v instead. 

• Phase transition model without source term. The physical solutions to the quadratic 
equation (|3.35 j) can be easily analyzed in the numerical study. Namely, we easily calculate 
that 

1 Dv St x(t i+ i) - x(tj_i) x(t i+1 ) - 2x(tj) + _ x(t i+1 ) - xjtj-t) 

v x Dt ~ 2 x(t i+1 ) - 2x(U) + x(ti-i) ' 5t 2 25t 

which is always positive. In light of the discussion at the end of Section we proceed 
as follows. If 

1 1 Dv 



2v x Dt 



v + - — 

is satisfied, then we compute q as 

-B(2Ap 3am -±%+-2v)+^A 



> VMam (4.49) 



pjam 

i o +- In tn ci n il i~\T n i I . , . , . 

v x Dt ^Hjam, 



(4.50) 



where A is given by (|3.36p . Thus, regardless what is the sign of v + — Apj am , the 
solution of p is always physical. 



Finally, the estimated quantities p, q depends on not only the types of traffic flow model and as- 
sumptions employed, but also the values of N. In order to compare different estimations on the 
same ground, we interpolate the continuous-time quantities such as v(t), w (i), p(t, x(t))and 
q(t, x(t)) with piecewise affine functions at the discrete points. In order to assess the estima- 
tion error, we will compute the relative L 1 error between the quantities obtained with N > 1, 
and the ones obtained with N = 1, i.e. the ground truth data. 

In the rest of this section, we will present the estimation errors associated with various 
traffic quantities including velocity and acceleration (Section 14.1 P ; p and q (Section 14.2)1 ; and 
the power demand function which is closely related to vehicle emission and fuel consumption 
(Section I4.3p . We will vary the value of N, or equivalently the sampling period of the vehicle 
trajectory, to see the impact on the estimation of first-order and second-order traffic quantities. 

Our study period is a fifteen minute interval from 4:00 pm to 4:15 pm, on a segment of 
1-80. A total number of 1724 vehicle trajectories are utilized. 



4.1 Reconstructing velocity and acceleration profiles 

In the first step of our study, we reconstruct the most basic traffic quantities: velocity and 
acceleration of the same vehcle, from (|4.41j) and (|4.42j) . Following the steps explained in 
the previous subsection, we summarize the mean and standard deviation of the relative error 
associated with velocity and acceleration below (Table 1). vt rU e refers to the ground-truth 
data, while vst is the approximated velocity by sampling every St seconds. 

From Table 1 we find that the estimation of velocity, which is a first-order quantity is 
much more accurate than the estimation of acceleration, which is a second-order quantity. 
Such a result can be intuitively confirmed by Figure [lj which shows the reconstructed velocity 
and acceleration trajectories by taking different sampling period. Since the time-trajectory of 
acceleration has worsened regularity compared to the trajectory of velocity, under-sampling 
(i.e. 5 = 3 seconds) tends to deteriorate the approximation quality of acceleration. But the 
velocity remains well-approximated in the L 1 sense. 
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Mean (%) 


Standard deviation (%) 


Sampling period St (s) 


0.5 


1 


2 


3 


0.5 


1 


2 


3 




Vtrue — Vst 


L 1 


3.27 


4.72 


6.36 


7.58 


1.50 


1.81 


2.31 


2.87 




Vtrue 






&true — 0>St 


L 1 


88.12 


98.77 


102.23 


102.53 


9.52 


3.10 


2.04 


2.04 




Q'true 


r. 1 



Table 1: Relative error using different sampling frequency. 
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Figure 1: Ground-true and estimated trajectories of velocity (left) and acceleration (right). 
4.2 Reconstructing p and q 

In this subsection, we will estimate p = pj a m — P, which is a first-order quantity, and q = 
q — q* , which is a second-order quantity. pj am and q* are fixed constants. A summary of the 
utilized dataset shows that the average speed of the 1724 cars under study is 17.5 mile/hour, 
which indicates that the traffic was in congested phase most of the time. Thus we employ 
three versions of the phase transition model with the congested phase: PTM with source and 
strongly stable traffic; PTM with source and less stable traffic; and PTM without source. 
In the velocity function (I3.14p . we choose A = 25, B = 1, pj a m = 1 (vehicle/m). The 
numerical results reported below are based on these values of parameters. We have repeated 
the numerical experiment with other values of A, B and pj am - The results show no qualitative 
difference. 



10 



4.2.1 PTM with source and strongly stable traffic 

The results of relative error are summarized in Table 2, with obvious meaning of notations. 
Figure [2] and [3] show the reconstructed trajectories associated with a given car. 





Mean (%) 


Standard deviation (%) 


Sampling period 5t (s) 


0.5 


1 


2 


3 


0.5 


1 


2 


3 




Ptrue — PSt 


L 1 


9.16 


10.65 


12.11 


13.02 


6.41 


6.38 


6.59 


6.80 




Ptrue 


J, 1 




Qtrue — Q8t 


L 1 


95.47 


100.36 


102.03 


102.08 


16.76 


5.64 


2.77 


2.56 




Qtrue 


L 1 



Table 2: Estimation based on phase transition model with source and strongly stable traffic. 
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Figure 2: PTM with source and strongly stable traffic: reconstructing the trajectory of p using 
different values of N. 
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Figure 3: PTM with source and strongly stable traffic: reconstructing the trajectory of q using 
different values of N. 
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4.2.2 PTM with source and less stable traffic 

The results of relative error are summarized in Table 3. Figure [Hand [5] show the reconstructed 
trajectories associated with a given car. 





Mean (%) 


Standard deviation (%) 


Sampling period St (s) 


0.5 


1 


2 


3 


0.5 


1 


2 


3 




Ptrue — PSt 


L 1 


12.17 


10.17 


11.13 


12.75 


16.36 


11.27 


10.28 


10.55 




Ptrue 


J A 




Qtrue — QSt 


L 1 


1161 


900 


772 


702 


1268 


977 


1110 


1039 




Qtrue 


r. 1 



Table 3: Estimation based on phase transition model with source and less stable traffic. 
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Figure 4: PTM with source and less stable traffic: reconstructing the trajectory of p using 
different values of N. 
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Figure 5: PTM with source and less stable traffic: reconstructing the trajectory of q using 
different values of N. 



We notice that for the phase transition model with source and less stable traffic, the error 
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for second-order quantity q is huge. This mainly steps from the fact that the L norm of the 
ground truth qtrue is very small. Recall (|4.47|) 



1 A 

1 A Dt 



3B v + 



lDv 
3 Dt 



and calculation (|4.48jl . We see that the above q is close (if not identical) to zero, even if slightly 
different approximations of v, v x and ^| are adapted. Such a fact contributes to the small L 1 
norm of qtrue- An example is shown in Figure [5j where the q almost vanishes everywhere. 



4.2.3 PTM without source 

In this model, the traffic estimates v, v x and j£ must satisfy condition (|4.49j) . We selected 
from the 1724 car trajectories the ones that satisfy such a condition. The number of such 
trajectories is 177. The results of relative error are summarized in Table 4. Figure [6] and [7] 
show the reconstructed trajectories associated with a given car. 





Mean (%) 


Standard deviation (%) 


Sampling period St (s) 


0.5 


1 


2 


3 


0.5 


1 


2 


3 




Ptrue — Pit 


L 1 


0.39 


0.56 


0.7 


0.75 


0.22 


0.28 


0.36 


0.42 




Ptrue 


J A 




Qtrue — QSt 


L 1 


5.41 


7.94 


9.87 


10.49 


2.79 


3.55 


4.37 


4.97 




Qtrue 


TA 



Table 4: Estimation based on phase transition model without source. 



From Table 4, we see that the estimations of both p and q are quite accurate. This is 
due to the fact that the velocity variations for these 177 vehicles are relatively small, and the 
traffic states are relatively stable around these vehicles. This can be partly seen from Figure 
El where the traffic density is quite stable along the vehicle trajectory. 




Figure 6: PTM without source. Curve fitting of p using different values of N. 
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Figure 7: PTM without source. Curve fitting of q using different values of N. 



4.3 Reconstructing emission and fuel consumption rates 



We use the Modal model (|3.38j) with 6 = to calculate the total tractive power requirements. 
Such a qua n tity i s closely related to vehicle emission and fuel consumption, as argued in 
Barth et al. ( 20071 ). The computation involves car speed v and acceleration 





Mean (%) 


Standard deviation (%) 


Sampling period St (s) 


0.5 


1 


2 


3 


0.5 


1 


2 


3 


Ptrue P&t j\ 


87.40 


97.60 


100.96 


101.26 


10.09 


3.29 


1.71 


1.63 


Ptrue jjl 



Table 5: Estimation of power demand. 

From Table 5, we conclude that due to the inaccurate approximation of acceleration, the 
estimation of engine power demand missed the target by almost 100%. 



5 Conclusion and Future Research 

In this paper, we showed the following. 

• Vehicle trajectory data (such as the ones collected using GPS) can be used to estimate 
Eulerian and Lagrangian traffic quantities such as v, j¥ and v x . The latter can fit into 
a second-order fluid model to estimate certain first- and second-order traffic quantities. 

• Overall, the first-order quantities such as velocity, density can be accurately estimated 
even if the sampling period of vehicle location is large (3 seconds). However, the esti- 
mation of second-order quantities such as momentum, acceleration, emission and fuel 
consumption rates deteriorates with long sampling period. In particular, GPS data alone 
is insufficient to estimate these quantities. 

• The phase transition model (PTM) with source and strongly stable traffic seem to be 
appropriate to treat various traffic scenarios. The PTM with source and less strongly 
stable traffic is not appropriate to estimate q. The PTM without source term seems to 
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approximate both p and q well, given that the traffic state is stable. In addition, it relies 
on quite restrictive assumption (|3.37p in order to deliver physical solutions. 

Some comments and suggestions for future study are as follows. 

• Colombo's original two phases approach fits traffic data well, yet the Siebel -Mauser 



source term allows the model to capture more complicate highways phenomenon (jSiebel and Mauser 
2006bh . 



The approximation of v x presented in this paper is more crude than those of v and Sf, 
and such approximation may lead to large error when the traffic is highly unstable. 

The data fitting depends on used measurements. If one uses only v and Dv/Dt is forced 
to assume strong local stability of traffic, thus neglecting some dynamics. On the other 
side the use of v x permits a better modelling but with rough data. 

The choice of traffic models depends on characteristics of the road segment under consid- 
eration (e.g. highways, arterial road and so on), the data quality (e.g. high/low sampling 
frequency), and the computational constraints. Intuitively, for long and stable highways, 
the phase transition model with source is a proper choice. For big urban arterial roads, 
phase transition model without source would be a good model. For short and unstable 
urban streets, one should use the classical LWR model. The next step is to test these 
models under various road conditions, as well as data and computational constraints. 
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